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In the most recent measurements of the reaction e+e“ —pp by the BABAR collaboration, new 
structures have been found with unknown origin. We examine a possible relation of the most distinct 
peak to the recently observed (j){2170). Alternatively, we analyse possible explanations due to the 
nucleon A and A A thresholds. The latter could explain a periodicity found in the data. 


I. INTRODUCTION 


The creation of nucleon-antinucleon {NN) pairs from the electromagnetic (em) current is one of the most 
fundamental baryonic processes and therefore an ideal ground to examine their interactions. This process 
exhibits the transition between the perturbative and nonperturbative regime of QCD. Close to the produc¬ 
tion threshold, it is dominated by meson exchange that gives rise to a strong enhancement Uj]. Increasing 
the center-of-mass energy, one finds an intermediate region, where this first thr'eshold enhancement has de¬ 
creased, though a description in terms of perturbative QCD is not yet possible. We will focus in particular 
on the region below the occurrence of the J/'h. This region is very poorly understood so far, so that a 
mainly conceptual analysis as performed here is appropriate. 

The nucleon form factors (NFFs) are the functions that parametrize the 'yNN vertex generated by the strong 
interaction. In the intermediate region of interest here, an effective proton FF shows several structures. We 
discuss for the first time the possible main sources of these structures. On the one hand, we examine whether 
the newly PDG-listed vector meson i;f)(2170) might be relevant here. On the other hand, we investigate if 
such an effect could be due to the final state interaction (FSI) at the thr'esholds of the first resonance excita- 
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tions, the NA + c.c. and AA thresholds, respectively. A possible AA-cu™ effect has been suggested by 
Rosner yD, in analogy to the case of pion photoproduction off the nucleon yl], but never examined further. 
For the treatment of the 0(2170), we consider a simultaneous description of the processes e'''e“ —)• NJV 
and eA eA, including observables for the proton and the neutron. Enhancements at the AA + c.c. and 
AA thresholds are treated in a simplified model calculation. 

The unphysical region of the NFFs is equally of interest, since a AA bound state, often denoted as baryo- 
nium, would be manifest below the physical threshold. Indications for such a state around an invariant mass 
of 1835 GeV have been found in the invariant mass-spectra of the decays J/'F —^ xpp, ^'(3686) —)• xpp 
with (x = 7 , uj, p, TT, p) and B'^ —)■ K~^pp, see e.g. Refs. IdJ-lTI]. The binding of a baryonium state could be 
generated by the final state interaction, giving rise to a pole below threshold, that could be accommodated 
in a recent FSI analysis The analytic continuation of the NFFs that we obtain from fits to data into the 
unphysical region is also examined with the help of logarithmic dispersion relations, including the region of 
a possible baryonium pole. Such an analytic continuation, however, requires the separation of the electric 
and magnetic NFFs over all the included kinematical range. This separation depends on the ratio between 
electric and magnetic FF. A higher precision for this ratio than from previous measurements is expected 
from the planned pp-annihilation experiment PANDA at FAIR IgL IlOll . 


For completeness, it may be worthwhile mentioning the findings of a partial wave analysis (PWA) of pp 


f. H 


elastic scattering in Ref. Illlll . where a peak in the ^D 2 partial wave has been explained as due to an S- 
matrix pole at 2.144 GeV and related to an unstable A A bound state. Also a recent PWA by the SAID 
collaboration lll 2 ll includes such pp data and finds a clear signal in the ^D 2 amplitude around this energy. 


Similarly, just below the AA threshold, a recent PWA including new neutron-proton scattering data il3h 
confirmed a pole related to the d*(2380), found at COSY ll^ . suggested to be a dibaryon. This raises the 
question whether similar mechanisms are at play in the pp system. 

The paper is structured as follows. In the rest of this introductory section, we give some basic definitions, 
discuss the contributions to the NFFs and experimental input. In Sec. HIl we examine a possible mani¬ 
festation of the 0(2170) in different nucleon observables. The inclusion of threshold cusps is examined 
in Sec. |III1 In Sec. |IVl we consider the analytical continuation of the NFFs to the unphysical region. We 
conclude with a discussion in Sec.lVl 


A. Definitions and prereqnisites 


For the description of the em process e~^{pi)e {p 2 ) —)• p{p 3 )p{p 4 ) we choose the center-of-mass (CM) 
frame, i.e. pip = {E,±ke) and P 3 4 = {E,±kp). The photon momentum q then determines the center- 
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of-mass energy by = {pi + p 2 )^ = ^CM ~ (2^')^- In our metric timelike q implies positive q"^. The 
three-momenta ke,kp appear in the phase-space factor /3 = kp/ke, which in the limit of neglecting the 
electron mass yields /3 « kp/E = — Arn^jq^, the velocity of the proton, and rup is the proton mass. 

We denote the emission angle of the proton by 6. The differential cross section in the one-photon-exchange 
approximation in this notation is 


da 

dQ 


0^(3 

4g2 




4m^ 


2m2 


(1-hcos 9)\GM{q )| H-^ sin e\GE{q )| 


( 1 ) 


where Ge and Gm denote the electric and magnetic Sachs form factors, respectively, and a = e^/(47r) = 
1/137.06 the fine-structure constant. C{q^) is the Sommerfeld-Gamow factor that accounts for the Coulomb 
interaction between the final-state particles 


G{q^) = 


TrarUr. 


i-e-y 


y = 


kr) 


Integrating over the full angular distribution gives the total cross section 


^e+e~^ppiq ) 


47ra^/3 


G{q^) 


3g2 

_ 47ra^/3 , 2 


2mi 


\GM{q^r + ^\GE{q^)\^ 


q^ 


3g2 


C{g‘) 11 + ^1 )l^. 


Thus, eliminating the kinematical factors from cr defines fhe effective form facfor Gefr 


( 2 ) 


( 3 ) 


|Geffl = 




\Gm\" 


\ 


1 + 


2m‘i 


( 4 ) 


For neufrons, fhe formulas are equivalenf excepf for fhe Sommerfeld-Gamow facfor which is nof presenf in 
fhaf case. Beyond fhe Coulomb FSI, higher order QED correcfions will be neglecfed in fhis work. The weak 
neufral currenf confribufion fo fhe measured cross secfion is also neglecfed. For fhe fime-reversed process, 
fhe phase space facfor is inverfed, yielding a^+^-^pp = /3^fTpp_j.e+e-■ 

Taking info accounf fhe angular dependence of pp producfion, one can express fhe differenfial cross secfion 
via fhe angular asymmefry A, 


wifh 


da 

dVL 


da 

dQ 


[1 -h ^cos^ 9], 

0 = 90 ° 


A 



-R^ 
+ i ?2 


( 5 ) 


( 6 ) 
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and determine from this explicitly the FF ratio R = \Ge/Gm\■ 

For many aspects, it is instructive to consider the vertex 7 pp in the helicity basis, i.e. the helicity-conserving 
Dirac and helicity-changing Pauli form factors Fi and F 2 , in order: 

^2 




4m| 


=i"l (<?') + 


On the one hand, this basis allows us to see directly that the threshold relation is by definition 


(V) 


GE(4mp) = GM(4mp). 


( 8 ) 


In addition, the asymptotic -dependence can also be conveniently given. For large —q^ = > 0, the 


Dirac and Pauli form factors can be predicted from perturbative QCD |[l5 


16h to behave like 


hm F,(Q2) = (g2)-(*+i) 
Q^^oo 


In 


A2 

/'-QCD, 


i = 1,2 


( 9 ) 


with 


7 = 2 - 1 -^ and /3 = 11 — ‘^Nr. 
' 3/ ^ 3 ^ 


( 10 ) 


Here, / is the QCD /3-function to one loop for the number of flavors Nf. The anomalous dimension 
7 depends weakly on the latter, 7 = 2.148,2.160,2.173 for Nf = 3,4,5, respectively. The analytic 
continuation of the logarithm to timelike momenta yields an additional term /A?) = \D.{q^/A?) — ivr, 

for g 2 > a 2 . 


For the asymptotic behavior of the form factors, we consider the Phragmen-Lindeloef theorem jlTII : “Let 


f{z) be an analytic function of z, regular and bounded in Im z > 0. If f{z) tends to the limits Li and L 2 
along the rays z = x + iO as x ^ ± 00 , then we must have Li = L 2 .” In particular, from this it follows that 
the imaginary part has to vanish in the asymptotic limit. 


Based on the reasonings o 


M 


' perturbative QCD including the analytic continuation into the timelike region 


1911 . a recent analysis 1201] of the proton effective FF includes a fit of the form 

A 


l^eff(9 )l — 


((/2)2(ln2(q2/A2) -F 7r2) ’ 


( 11 ) 


with the parameters from a fit to data prior to a recent measurement by the BABAR collaboration 12111 . 
given as A = 72 GeV“^ and A = 0.52 GeV. 
















5 



a) b) c) 


FIG. 1: The different types of diagrams contributing to the NFFs above the pp threshold, as discussed in Sec. IIBI The 
wiggly line denotes the photon, the thin solid line the (anti)nucleon, the thick solid line an (anti)nucleon resonance 
and the dashed line denotes all possible mesons, e.g. pions in a)+b). 

B. Possible contributions 

In the energy-regime below the J/'h, we will focus on possible hadronic contributions to the NFFs. The 
size of the different possible contributions is unknown. Therefore we consider here the individual diagrams 
in a pioneering study, neglecting interference effects between them. 

Three Feynman graphs representing the different types of diagrams contributing to the NFFs are shown in 
Fig. [U The first two refer to a mainly baryonic, the last to a mesonic contribution. Diagram a) represents 
the final-state interaction (FSI) in form of meson exchange. The pion exchange shown can be replaced by 
any number of suitable mesons. Diagram b) shows one possible excitation of a resonance, e.g. a A, in the 
FSI diagram. Possible re-excitations are assumed. Diagrams of type c) are usually ignored above the pp 
threshold. However, in general they can contribute, e.g. from the (/>(2170). 


C. Data 

Experimental information on the proton FFs for timelike momenta is available from a number of measure- 




FIG. 2: The processes mainly considered for NFF determinations. Left: Annihilation measured at LEAR, planned at 
FAIR 1^ 3 Right: Production with ISR, measured e.g. at BABAR. 


ments dating back to 1976, for a detailed discussion, see Ref. i23h . However, the only ones that inc 


sufficient precision on the angular distribution to disentangle Ge and Gm, are those from LEAR 12511 . 


ude 


BABAR i2in . see Eig. [2l and BESIII 12611 . In fact, the first dispersion analyses of timelike nucleon form 
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FIG. 3: The total cross section a(e~^e —?■ pp) measured at BABAR, depending on the invariant mass Mpp. The 
initial state radiation is factorized out and the J/'I' and T'(2s) peaks have been removed 1^ . 



H], 


FIG. 4; The effective form factor of the proton measured in the process a{e^e PPl) at BABAR 
shown with different fits from Ref. l20l] . given by a conventional dipole-monopole form or the pQCD parametrization 
Eq. O- 


factor data 1241] were hampered by a missing separation oi Ge and Gm- The recent BESIII measurement 


yields good agreement with BABAR, in particular confirming an enhancement at 2.4 GeV. However, the 
emphasis of BESIII lies above the region of our main interest here. Due to the much higher precision of the 
BABAR data compared to the ones from EEAR, in this work we focus on the former, for the corresponding 
cross sections see Eig. [3l The ridge and bump structure occumng here are also inherited by the effective 
proton EE that is shown in Eig. |4l The insufficiency of a pQCD description is shown by the fit Eq. (fTTl) from 
Ref. 1^ . where it is argued that the oscillations around the pQCD fit are due to ESI. 
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In this work, the angular information is included in the form of the FF ratio \Ge/Gm\ as provided by 
BABAR Q. 

The neutron FFs for timelike momenta have been measured in the process e'''e“ —)• nn by the FENICE 


collaboration 12711 and recently at the VEPP-2000 collider with the “spherical neutral detector” (SND) 12811 . 


As for the proton case, the precision of the latest neutron measurement exceeds by far that of the preceding 
experiments. However, the accuracy does not suffice to determine the neutron form factor ratio. 


In the spacelike region, we include explicitly the mostprecise differential ep scattering cross sections 12911 
and in addition the polarization world data (see Ref. OOll for a list). Eor the neutron, we want to emphasize 
that the phenomenological extraction of neutron FFs yields larger uncertainties than for the proton case, 
since there is no free target. The measurements require light nuclei targets like or ^He for quasi-elastic 
scattering, for details see e.g. Ref. 00(1 . Here, we use the data on and G^ for a better visibility 
compared to the equivalent use of G^ and the ratio. 


II. THE PROCESS e+e” ^ NN AND THE <()(2170) 


In this section we consider the possible contribution to the NFFs from the 0(2170), corresponding to di¬ 
agram c) in Fig. [T] This refers to a structure that has been found in different processes and that is at 
the moment classified by the PDG as the only light unflavored vector meson above the A^iV-threshold. 
In the following, we use the PDG notation even though some of these structures have also been denoted 
as X/Y (2175) or simply as cj)”. As shown in Tab. IH measurements at BES, BABAR and BEEEE have 
found signals in this mass region, albeit over some interval. Even the central values spread in the range 
2.08 — 2.19 GeV for the mass and 58 — 192 MeV for the width. This might correspond to the uncertainty 
of the separation from non-resonant background and/or the possible existence of multiple interfering reso¬ 
nances in this range. Also the isospin is given as definite, yielding altogether = 0“1 . However, 

the assignment of quantum numbers should be taken with a grain of salt, mainly due to the limited statistics, 
for a discussion see for example Ref. H. Different suggestions about the origin of the 0(2170) have been 


put forward. It has been interpreted, for example, as a tetraquark state ll32ll . a hybrid ssg resonance 1133h 
or to a large extent as a 0(lO2O)iF.A state 0411 . It can also be generated in a chiral Eagrangian approach 


for 0(lO2O)//o(98O) S-wave scattering by their self-interactions OSh . However, to our knowledge it has 
not been considered in relation to the NEEs. We will do so, first by focusing on the effective proton EE and 
second in a simultaneous treatment of different proton and neutron measurements for space- and timelike 
momenta. The individual form factors that correspond to a Breit-Wigner resonance structure with mass 
and width F^ behave like F{q^) oc 1/(M| — so that the effective FF can be fitted to the 











absolute value of the latter. 


process 

mass (MeV) 

width (MeV) 

J/T- ^ r](j)fo{980) [BES] 

2186 ± 10 ±6 

65 ± 23 ± 17 

—>■ [BABAR] 

2125 ± 22 ± 10 

61 ± 50 ± 13 

e+e" ^ K+K-nnj [BABAR] 

2175 ± 10 ± 15 

58 ± 16 ± 20 

e+e" A+A-TT+TT-y [BELLE] 

2079±13t^^ 

192 ± 23^61 

e+e" ^ A+A-TT+TT-y [BABAR] 

2192 ± 14 

71 ±21 

e+e" ^ A+A-TrOrrOy [BABAR] 

2169 ±20 

102 ± 27 


TABLE I; Observations of the ^(2170) from the review of particle properties Obll . 


A. Individual fits 

A description of the proton effective form factor has been attempted by several groups over the years with the 
main emphasis on either the perturbative QCD part or the impact of vector mesons below the pp threshold. 
Neither these nor conventional fits of dipoles, monopoles or products of these can fully accommodate the 
structures in the currently most precise relevant cross sections or in see Figs. [3j IH In the overall 




FIG. 5: Illustrative pre-fits for comparison to PDG-given (^(2170), using only GeS-data. Including a resonance at 
Ml = 2.125 GeV with a width of Ei = 0.09 GeV (left) and Ei = 0.33 GeV (right). The fit on the right also includes 
a resonance at M 2 = 2.43 GeV and E 2 = 0.23 GeV. 

concave function the first structure occurs as a mainly convex part for invariant masses of around 
2 — 2.25 GeV. A satisfactory fit to the data requires to take this structure into account in some way in the 
parametrization. However, this is a prime example of the ambiguities that can appear in separating a possible 
resonance structure from the background. For comparison with PDG-values, see Tab. HI we perform test fits 
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to Ggg with terms that correspond to Breit-Wigner (BW) shapes in the cross section and 5 effective pole 
terms below threshold for the background description. We find that a large range of values of masses and 
widths in this region strongly improves the data description, M ^ (2—2.18) GeV and F Ri (0.05—0.5) GeV. 
In Fig. [5l we show examplary fits with a narrow resonance on the left and a broader one on the right. 
Remarkable is also a structure peaked around « 5.9 GeV^. Allowing a second resonance in this region 
yields one at M = 2.43 GeV with a width of F = 0.23 GeV. This cannot be regarded as a rigorous analysis, 
since these resonance structures largely overlap and are not separable from the background either. However, 
it is undebatable that an additional structure peaked around M ^ 2.43 GeV improves the data description. 


B. Simultaneous fits 


In this section, we combine the fits from the last section with more available data on the NFFs, see 
Sec. lICI for references. These data comprise 7 different sets, 4 for the proton and 3 for the neutron. For the 
proton, we consider the differential cross sections and the ratio Ge/Gm from polarization observables on 
the scattering side in addition to the effective FF and \Ge/Gm\ on the production side. For the neutron, 
we include Ge and Gm from scattering data and again the effective FF on the production side. In order 
to weight the different data sets equally, their impact on the -function to be minimized is determined by 
their number of data points. Still, to avoid a dominance by the other sets, we fix the mass of the resonance 
structure at M = 2.125 GeV. The width obtained in the simultaneous fit is F = 0.088 GeV. The larger 
number of data points compared to the previous section requires more effective pole terms in the unphysical 
region, in particular for a separation of the isospin channels due to the inclusion of the neutron. As the basic 
framework, we proceed in a similar way to our analysis of spacelike NFFs 1^ . To be specific, we include 
parametrizations of the 27r, KK and pvr-continuum, the lo- and (/(-contribution and effective pole terms. 
Here, the latter are restricted to the region 1 GeV^ < < 3.52 GeV^ and limited to a number of three in 

the isoscalar and five in the isovector channel. 

We show the fit results in Figs.[6l|7l[ll|9l Note that the first thr'ee subgraphs refer to the timelike and the 
remaining four to the spacelike region. For the proton effective FF, we included a second set of BABAR 
data at higher energies, that are still well described 


ry our fit. For the neutron we included the recent 
measurement from the VEPP-2000 e'''e“ collider 112811 . In the range from threshold up to = 4 GeV^, 


their -dependence is very similar to the proton case. Only above this, two further data points given by 
the FENICE collaboration indicate a less steep fall-off. Unfortunately, the neutron data is too sparse in this 
region to constrain a possible manifestation of the (/)(2170) around = (2.125 GeV)^ « 4.5 GeV^. Also, 
the current level of statistics of the angular distribution of nn events is too low to determine the | G^ / G^ \ 
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FIG. 6: The effective FF of the proton (left) and the neutron (right). G^g is complemented here by the subsequently 


published higher-energy data 


published data from FENICE JZ 


3811 . The recent SND measurement 12811 is in good agreement with the only previously 



q'^lGeV'^] 



EIG. 7: The form factor ratio of the proton for space- (left) and timelike (right) momenta from the combined fit to 
space- and timelike data. 


ratio ll28ll . 


So we are left with the ratio to search for direct indications of a resonance at ^ 4.5 GeV^. 

Indeed, a slight dip occurs here, see Fig. |7J as soon as we include the resonance term. The same quantity in 
the spacelike region is also well reproduced. For the neutron FFs in the spacelike region, see Fig. [8l we want 
to emphasize the sizeable uncertainties in their extraction from electron scattering off light nuclear targets 


like or ^He, for details see e.g. Ref 


. H. Fi 


Finally, the electron-proton scattering cross sections shown 


in Fig. |9] are by far the largest data set with 1422 out of the total 1627 points. Therefore weighting each set 
equally disfavors the larger sets in a sense, still giving reasonable agreement in this case. One could also 
here refit the normalization individually for different parts and treat the uncertainties as discussed in detail 
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FIG. 8: The neutron electric (left) and magnetic (right) form factor from the combined fit to space- and timelike data. 


in Ref. 



However, for a more conceptual work like the one at hand we refrain from such a procedure. 


HI. THRESHOLD ENHANCEMENTS 


In this part, we outline an alternative origin of the structures found in the data by carrying out fits to 
Ggg only. Remarkable in the last section is the position of the peaks, or rather kinks, that are necessary 
to improve the fits. The positions can be chosen as the threshold energies of the pA + c.c. and the AA 
states, or pp2TT and pp^Tr, respectively. The occuring A resonance would emit a pion (> 99%), or a 
photon (< 1%). Since the backgrounds to the pp final states are subtracted, as discussed in detail by the 
BABAR collaboration 1^ . the pion or photon has to be absorbed by the other (anti)baryon. In general, the 
interaction between the final states can of course comprise further exchange of on- or off-shell mesons which 
should be treated systematically in an effective field theory (EFT) framework. Close to the pp threshold, 
the final-state interaction (ESI) can be computed via a Eippmann-Schwinger equation. Such a procedure 
has recently been updated in Ref. Uj], employing a static interaction ^tential derived in chiral EFT 091] 


that has been fitted to a partial wave analysis of pp scattering data idOT . Moreover, this is based on the 


assumption of a real and constant bare vertex function. In the region of validity of such an approach 
{Mpp — 2mp) < 0.1 GeV, the decrease in can be well reproduced. Beyond this range, a calculation 
in this framework breaks down. However, one naively expects the interaction to further decrease due to the 
increasing relative velocity of the two final states. In such a naive reasoning, the excitation of a resonance 
would lead to the same threshold kinematics, just shifted in energy, and thus could give rise to an enhanced 
ESI. After the decay of the A resonance, which dominates here, one would be left with mainly the pp ESI. 
For the possible necessiW to resum the loops corresponding to re-excitations, future work could proceed 


along the lines of Ref. 14ill . This might allow us to distinguish the possible origins of the structures found in 
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0 20 40 60 80 100 120 140 

scattering angle 9 [°\ 


FIG. 9: The electron-proton scattering cross sections from MAMI ll29ll . On the left, the energy of the incoming 


electron is given, which together with the scattering angle determines Q^. 


Ggjj. For the moment, however, we merely illustrate some possible contributions from the triangle diagram 
shown in Fig. [TO] Also these are only roughly approximated due to the lack of information on the vertices, 
in particular their momentum dependence. 


A. Inclusion of the NA + c.c. and A A thresholds 


We consider the triangle graphs with virtual NAtt, see Fig.[T0l and AAtt in order to approximate possible 
cusp effects. However, the vertices are not well known for these kinematics. While, e.g., for the ANtt 
transition the coupling constant at = —M^ are known, the form factors and their general dependence 
on the different momenta is all but well known. This would be relevant if we were to evaluate the triangle 
diagram in full glory, which we do not attempt here. What can be obtained most easily though, is the scalar 
part of the integral. This is proportional to the analytically well-known Passaiino-Veltman integral C'o(k). 
As defined in Ref. |^, (slightly different conventions), and implemented in LoopTools |^], this depends 
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P3 


PA 


FIG. 10; Triangle graph with virtual IVAtt state. Notations are as in Fig.[T] 


on the configuration («;) of virtual particle masses and external particle four-momenta: 

If d'^k 

Co{K=Pk,Pl,mi,m2,m3) = ^ J _ ^ 2 ] ( 12 ) 

with Ki = P 3 ,Pa, Mtt, mAT, ttia and K 2 = Ps,Pa, MT^,mA,'iTLA for the two cases considered here, omitting 
the ie. The inclusion of the A width also changes the analytic structure. In the following, we denote the 
configurations corresponding to ki, K 2 with A widths by tci, ti; 2 - We show the absolute value and imaginary 
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FIG. 11; The q^-dependence from the scalar Passarino-Veltman triangle diagrams with virtual NAtt and AAtt states 
compared to from Ref. . 


part of the configurations ki and K 2 in Fig. [TT] Also shown is how the inclusion of the A width partly 
smears out the cusp effect. Taking the loop momenta in the numerator into account, one can reduce the 
graph to a sum of n-point functions with n < 3. The momenta only partly cancel against those in the poorly 
known form factors, so that an additional smearing of the result is expected. 


B. Cusp fits 

In this section, we show how the scalar parts of the relevant triangle diagram compares to the convex 
structures in Even after inclusion of the width, the remaining enhancements have the right position 
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and shapes to improve a pure pole fit. As in Sec. Ill A[ we fit only include 5 effective pole terms (below 
threshold) and now replace the explicit resonance terms by the loop structures from the last section. In order 
to account for the form factors at the vertices and a smearing as discussed in the last section, we include one 
form factor for each loop 


Hq") = 


(13) 


1 + 

with the respective fitted cut-off parameter. Additionally, the overall size of the loop contributions 

is allowed to vary by a fit parameter /ata/aa for o^oh loop. The fit result with /ata = 0.02, /aa = 0.3 



FIG. 12: Fit including effective pole terms and the scalar parts of the triangle diagrams with virtual NAn and AAtt 
states to Ggjj from Ref. 1^ . 

and Aata = 10 GeV, Aaa = 1-7 GeV is shown in Fig. [121 The fit parameters / na/aa of natural size, 
as expected. 


IV. THE UNPHYSICAL REGION 


In this section, we discuss the NFFs in the region of fg = 4M^ < t < tph = which is not accessible 
by direct measurements, but by analytic continuation in f = —Q^. An additional particle emission 
from the initial state _proton can lower the energy of the (virtual proton) to reach below the threshold, 
as discussed in Ref. jddl] for the process pp e~^e~nQ. Without model assumptions though, one can 


relate the information in this region to the physical ones by means of dispersion relations (DRs). In many 
applications of these, the higher energy parts of the spectral function are not particularly relevant or are 
suppressed by subtractions. However, it can also be of interest to use (experimental) information from the 
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physical timelike region which only determines the absolute value of the NFFs. To this aim, it is instructive 
to use a DR for the logarithm, see e.g. Refs. |3-49]. In principle, this also allows for a separation of the 


FF phase 6{t) and modulus in the representation G{t) = The relative phase of and Gm 

in turn, can be obtained in polarization measurements, as planned for PANDA at FAIR. This phase might 
help to understand the origin of the structures in Moreover, with ideally accurate data in the space- 
and timelike physical region one could obtain information on both the modulus and phase of the FFs in the 
unphysical region, including the latter above production threshold. 

One can start from a subtracted DR for the function ln[G(f)/G(0)]/(f-v/fo ~ f)- For f < 0, we evaluate the 
DR 


InG(f) = InG(O) + 


- t 


TT 


/■“ In I 

Jto F(F 


ln|G(t')/G'(0)| 


poo 

■.dt'= / , 

Jtn 


(14) 


'(F - t)y/t' - to 

where the first term vanishes due to the normalization G£;(0) = Gm{^)/ d-p = 1- Experimental information 
on this integral equation (fT4l) is available in the spacelike region f < 0 on G{t) and in the timelike region 
for t > tph on the modulus |G(f)|. One can thus split the integral into the known part above t > tph and 
the remaining part with unknown integrand, as suggested in Ref. 1471] . The resulting integral equation is 


commonly denoted as an inhomogeneous Fredholm equation of the first kind ll50n . In general, the solution 


for the unknown part of In |G(f)| can be searched for by discretizing the integral. Our first choice would be 
a number of discretization points equal to the number of input points from the physical region, giving a set 
of n linear equations with n variables. However, the problem is strongly ill-conditioned, with small changes 
in the input leading to large changes in the solution. Therefore, additional information is required to solve 

^ and consider the integral contributions 


the original integral equation. We proceed similar to Refs. 11471 


to the logarithm In |G(f)| in the spacelike region, using definite values for the known part above tp^ 

C^ph 

I tph ^to 


p'oo piph 

\nG{t) — / I{t,to,t')dt' = / I{t,to,t')dt', t < 0. 

j t^h J tc\ 


(15) 


In contrast to Refs. 1147 , 


4811 ■ we use as input for the Ihs of Eq. (fTSl) our discretized result of a simultaneous 


fit to data in all accessible regions, see Sec.ITTl and obtain as an example fit the result shown in Eig.[T3l We 
minimize the difference between left- and right-hand-side of Eq. (fTSl) . while also limiting the total curvature 
of the FFs in the unphysical region, r = {{d?‘\G{t)\/dt^))'^dt. For the result, we find a large dependence 

on fhe sfrengfh of fhe curvafure limifafion and on the range that we use for input from f < 0. Thus reliable 
errors can not be quantified here. 

However, we wanf fo point out that large enhancements in the FF modulus just below production thi'esh- 
old are possible. In the example shown in Fig. [131 one can see such an enhancement just below tph, as 
expected in the case of a baryonium pole. For the planned precision of the forthcoming measurements at 
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FIG. 13; An exemplary result for the modulus of electric and magnetic form factor obtained from the logarithmic 
integral Eq. dl. The NFFs for t < 0 are considered via the differential cross sections, therefore the data cannot 
be shown directly in this form. For t > tph, \Gm\ is closer to |Geff| than as expected. We note that large 

enhancements in the FF modulus just below production threshold are possible and should be evaluated via this method 
with future PANDA data. 

PANDA, these possibilities should be further evaluated. Encouraging in this regard are also previous results 
for the pion FF 1^ . where the predictive power of this method can be impressively illustrated. A major 
source of complication in the nucleon case is the neccessity of two FFs and their separation. The emphasis 
on measuring the angular distribution at PANDA will have particular impact on this separation. 


V. DISCUSSION AND CONCLUSIONS 

In this paper, we have examined the em pp creation and the mechanisms that dominate it in the domain of 
nonperturbative QCD. Specifically, we have discussed possible contributions to the NFFs corresponding to 
a vector meson recently listed in the PDG as cj){2l70) or from FSI at the NA + c.c. and AA thresholds. 

We have included the (j){2l70) in simultaneous fits to proton and neutron FFs for space- and timelike mo¬ 
menta and found good agreement with the existing data. In particular, we included recent measurements 
on the neutron effective FF. In contrast to the previous FENICE experiment and analyses of this, the recent 
SND data shows a very similar behavior to the proton effective FF over a large range, which we can describe 
well in our approach. However, the range around the cj){2l70) calls for further neutron measurements to 
allow for a determination of the isospin channel of the structures in 

It may be worthwhile mentioning here, that similar fits to data only in the spacelike region, as performed 
in Ref. 1^ . found exactly two (“effective”) poles in the physical timelike region, one at 2.14 GeV and one 
at 2.4 GeV, each 10 — 20 MeV below the N A and A A thresholds, respectively. Accordingly, we have 
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also examined possible contributions from the final state interactions at these thresholds. Taking approx¬ 
imations for the FSI into account allows for a similarly good description of the data as the inclusion 
of the (/)(2170). The occurrence of peaks in around both A thresholds might favor this explanation. 
However, future calculations should include the singularity structures of any possibly contributing diagram 
and all interferences. In particular for the case that some structures indeed exist below the thresholds, a 
resummation of the FSI diagrams is clearly necessary to calculate the pole of such a bound state. In this 
context it may be of interest, that the small deviation at ~ 2.25 GeV lies close to the A A threshold. 

In order to distinguish between the possible effects we are awaiting polarization measurements at FAIR 
from which one can extract the relative phase of Ge and Gm- This will also improve the precision of the 
analytic continuation to the region of a possible baryonium pole. Using logarithmic dispersion relations, we 
have found that such a pole with large contributions to the NFFs could exist. 
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